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INTRODUCTION 


The  initial  report  on  the  MINIMUF-3  model  (ref  1)  found  widespread  interest 
throughout  the  hf  community.  That  model  has  been  successfully  adapted  for  use  on  a 
variety  of  minicomputers  and  microcomputers  and  has  been  used  in  numerous  ways.  As 
long  as  the  800-8000  km  range  constraint  discussed  in  ref  1 is  not  exceeded,  the  original 
MINIMUF-3  will  perform  within  confidence  limits  specified  in  the  first  verification  study. 
But  the  development  of  new  applications  involving  very  long  transmission  paths  uncovered 
several  bugs  that  had  to  be  corrected.  The  changes  made  do  not  affect  the  original  accuracy 
limits  but  do  allow  MINIMUF  to  be  used  out  to  the  antipodal  point.  The  purposes  of  this 
document  are  to  present  the  most  current  model,  called  MINIMUF-3. 5,  to  describe  the 
improvements  made,  and  to  present  a data  set  with  which  each  user  can  check  out  his  ver- 
sion of  MINIMUF-3. 5.  The  original  report,  ref  1,  explains  the  theory  behind  this  model. 


CHANGES  TO  ORIGINAL  MINIMUF  MODEL 

The  model  listed  in  the  BASIC  program  on  page  41  of  reference  1 contained  two 
problems,  as  follows. 

1.  Line  number  891  is  a conditional  branch: 

891  IF  COS(L0+Y2)>-0.26  THEN  900 

If  the  branch  to  line  900  is  not  taken,  subsequent  execution  of  the  line 

1 220  G2=G2*(  1-0. 1 *EXP((K9-24)/3)> 
references  the  variable  K9,  which  will  not  have  been  properly  defined. 

The  problem  is  eliminated  by  insertion  of  the  following  line  of  code: 

892  K9=0 

2.  For  very  long  path  lengths,  the  computed  MUF  is  unrealistically  low.  For 
paths  sufficiently  long,  a fatal  error  can  occur  at  line  898  or  line  1000  with  a negative  argu- 
ment for  the  square  root  function.  This  is  because  the  value  of  M9  (computed  in  line  897 
or  990)  reaches  a maximum  value  for  a path  length  of  about  8000  km  (the  longest  path 
used  in  the  verification  tests)  and  then  decreases  in  magnitude,  becoming  negative  for  path 
lengths  of  about  16  000  km.  The  following  changes  correct  this  problem: 

Delete  lines  897,  898,  990  and  1000.  Insert  the  lines 

897  M9=S1N(2.5*G1  *K5  MIN  PI/2) 

898  M9=1+2.5*M9*SQR(M9) 

990  M9=S1N(2.5*G1  *K5  MIN  PI/2) 

1000  M9=1+2.5*M9*SQR(M9) 

In  addition  to  these  changes,  MINIMUF-3. 5 uses  a somewhat  different  solution  for 
the  spherical  trigonometry  problems  involved  in  computing  path  length  and  control  point 
coordinates.  This  change  simplifies  the  code  and  eliminates  the  need  for  GOSUBs  to  exter- 
nal subroutines.  It  affects  neither  the  solution  obtained  nor  the  computed  value  of  MUF. 


*NOSC  TR  186.  MINIMUF-3:  A Simplified  Hf  MUF  Prediction  Algorithm,  by  RB  Rose,  JN  Martin,  and 
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REVISED  MINIMUF  BASIC  PROGRAM 

A listing  of  the  MINIMUF-3.5  program  is  included.  Lines  100  through  720  contain 
a small  driver  which  allows  the  model  to  be  exercised.  The  actual  MINIMUF  program  starts 
at  line  1000. 

The  input  variables  for  the  MINIMUF  program  are  as  follows: 

LI  - Transmitter  latitude,  radians  (-rr/2  < LI  < jr/2) 

W1  - Transmitter  west  longitude,  radians  (-27T  < W]  < 2it) 

L2  - Receiver  latitude,  radians  (-jt,/2  < L2  < n/2) 

W2  - Receiver  west  longitude,  radians  (-2rr  < W2  < 2ir) 

M0  - Month  (1  <M0<  12) 

D6  - Day  ( 1 < D6  < 3 1 ) 

T5  - Time  (UT).  hours  (0.0  ^ T5  < 24.0) 

J9  - Output  MUF.  MHz 
S9  - Sunspot  number 
PI  - 3.141593 
P0  - 1.570796 


In  using  the  test  driver  program,  transmitter  and  receiver  coordinates  should  be 
input  in  degrees;  the  driver  converts  these  to  radians  as  required  by  the  MINIMUF 
subroutine. 

Figure  1 is  a sample  output  listing  that  users  may  find  helpful  in  getting  their  version 


of  MINIMUF  to  work. 


TRAMSNITTER  LAT.  LON  • 21, 
RECEIVER  LAT,  LON  ■ 30,122 
DATE  <0AV,N0NTN>  - 17,10 
SUNSPOT  MUNAER  • 100 


IK 


DATE:  17  OCT 
TRAHSH1TTER  LOCATION: 
RECEIVER  LOCATION: 
SUNSPOT  NUNOER  ■ 110 


NOUR 

NUF(HHZ) 

0 

32.0 

1 

32.0 

2 

32.0 

3 

29.9 

4 

23.0 

3 

22.0 

t 

20.9 

7 

19.3 

0 

10.0 

9 

10.9 

10 

16.0 

11 

13.2 

12 

14.0 

13 

14.1 

14 

13.7 

13 

21.0 

10 

27.6 

17 

31.3 

10 

32.0 

19 

32.0 

20 

32.0 

21 

32.0 

22 

32.0 

23 

32.0 

PRESS  RETURN  TO  PERF 

LATITUDE  21.00, 
LATITUDE  38.ee, 


LONGITUDE  1SS.I 
LONGITUDE  122.1 


Figure  1 . Example  of  a 24-hour  MUF  lisiing  of  MINLMUF-3.5. 
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APPLICATION  TIPS 


With  increasing  solar  activity,  user  interest  in  updating  MINIMUF  to  reflect  current 
solar  activity  has  also  increased.  The  updating  method  found  to  be  most  effective  was  to 
vary  the  sunspot  number  input  parameter  as  a function  of  107  mm  solar  radio  flux.  Because 
of  the  lag  in  F-region  response  to  a rapid  increase  in  solar  activity,  it  is  best  to  use  either  a 
5-day,  1 5-day  or  90-day  running  mean  average  of  the  1 07  mm  flux.  The  type  of  application 
will  determine  which  is  best.  The  5-day  mean  is  a shorter-term  more  dynamic  input,  while 
the  90-day  mean  is  more  applicable  for  long-term  planning.  These  flux  values  can  be 
acquired  from 

Space  Environment  Services  Center,  Boulder,  Colorado 

WWV  transmissions  at  1 8 minutes  after  each  hour. 

The  conversion  from  107  mm  flux  to  sunspot  number  is  accomplished  by  the  graph  shown 
as  figure  2. 

Two  other  points  are  borne  out  by  field  testing.  First,  MINIMUF  is  an  F-region 
approximation.  Any  intervention  by  E-region  modes  of  propagation,  either  as  multiple  E 
or  EF  complex  modes,  is  not  predictable  on  MINIMUF.  Such  operational  situations,  how- 
ever, are  proving  to  represent  only  a small  percentage  of  the  total.  Second,  MINIMUF  has 
the  greatest  accuracy  within  the  one-  and  two-hop  ranges,  between  about  800  and  8000  km. 
Predictions  for  transmission  paths  longer  than  this  should  be  used  with  some  caution. 


Figure  2.  Relationship  between  smoothed  mean  Zurich  sunspot  number 
and  smoothed  mean  107  mm  solar  flux. 
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MINIMUF-3.5  PROGRAM 


1000  REN 

1010  K7-SIN<L!>tSINa2>+C0S<Ll)tC0SCL2>tC0S<M2-Hl> 

1020  G1-ACS<K7  NAX  -1M.0E-3  NIN  l-t.OE-8) 

1030  K6M.39tGl 
1040  K6«K6  NAX  1 
1050  K3*1'K6 
1060  J9M80 

1070  FOR  Kl-l'<2tK6>  TO  1-1/<2*K6>  STEP  0.9999-1/K6 
1060  IF  K5-1  THEN  1100 
1090  K5-0.5 

nee  p»sin<l2> 

1110  Q«C0S<L2> 

1120  A«<SINai>-PtCOS<Gl>>/(Q*SIN<Gl>> 

1130  B-GUK1 

1140  C«P*COS<B>+QtSIN(B>*A 

1150  D»<C0S<B>-C*P>/<Q*SQR<1-C*C>> 

1160  D«ACS(D  NAX  -l+l.OE-5  NIN  1-1.0E-8) 

1170  H0-M2*SGN<SIN<U1-H2))*D 
1180  IF  MO«>0  THEN  1200 
1190  M0*M0+P1 
1200  IF  H0<P1  THEN  1220 
1210  M0*M0-P1 

1220  LO«P0-ACS<C  NAX  -1M.0E-5  NIN  1-1.0E-8) 

1230  Yl*0.0172$< 10><H0-1)$30.4+D6> 

1240  Y2«0.409*CO$<Y1> 

1250  K8»3.82*H0+12^O.13*<SIN<Y1>M.2*SIN<2*Y1>> 

1260  K8-K8-12t<nSGNCK8-24>>*SGN<ABS(K8-24>) 

1270  IF  COS (L0+Y2) >-0.26  THEN  1350 

1280  K9-0 

1290  GO*0 

1300  H9-2.5SG10K5 

1310  N9-H9  NIN  PO 

1320  H9*SIN<N9> 

1330  H9-l42.5tN9tSQR<N9> 

1340  GO  TO  1590 

1350  K9-(-0.264SIN(Y2>*SIN(L0>>/<C0S(Y2>SC0S(Le>-n.eE-3> 

1 360  K9» 1 2-ATN ( K9/SQR  < ABS  < 1 -K9*K9 > > > *7 . 639437 
1370  T »X0-K9^2^ 1 2 1 ( 1 -SGN < K8-K9/2 > >*SGN< ABS<K8-K9/2>  > 

1 380  T4«K0^K9/2- 1 2* <1 ♦SGN  < K8+K9/2-24 > > $ SGN < ABS  < K8+K9/2-24  > > 
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1398  C8~ABS(COS<L0+Y2>) 

1400  T9*9.7<C0t9.6 
1410  IF  T9>0. 1 THEN  1430 
1420  T9*0. 1 
1430  H9S2.5<G1<K5 
1440  H9*H9  MIN  P0 
1450  H9«SIN<H9> 

1460  H9«1*2.5<H9<SQR<H9> 

1470  IF  T4<T  THEN  1500 

1480  IF  <T5-T>«T4-T5>>0  THEN  1510 

1490  GO  TO  1640 

1500  IF  <T5-T4>XT-T5>>0  THEN  1640 

1510  T6*T5*H2X1*SGN<T-T5>XSGN<ABS<T-T5>> 

1520  G9*PIXT6-T>/'K9 
1530  G8<PKT9/K9 
1540  U*<T-T6>/T9 

1550  G0*C0XSIN<G9>+G8XEXP<U>-COS<G9>>>/<1+G8<G8> 

1560  G7<C0XG8«EXP<-K9/T9>+  1 > ><EXP< <K9-24)/2>/< UG8<G8> 
1570  IF  G0* >G7  THEN  1590 
1580  G0*C7 

1 590  G2«  < 1 *S9s250  > <N9<SQR  < 6+58<SQR  < GO  ) ) 

1 60O  G2*G2< <1-0. 1<EXP< <K9-24>/3> ) 

1618  G2*G2« 1 ♦< 1 -SGN<L1 ><SGN<L2>  ><0. 1 > 

1620  G2«G2«l-0.  1<<1+SGN<ABS<SIN<L0))-COS<L0))>> 

1638  GO  TO  1700 

1 640  T6*T5+ 1 2X 1 ♦SGN  < T4-T5  > > <$GN < ABS < T4-T5 ) > 

1650  G8-PUT9/K9 
1660  U«<T4-T6>/2 
1670  U1*-K9'T9 

1680  G0<C0V<G8XEXP<U1  )+l ) )<EXP<U>s<  1+G8<G8> 

1690  GO  TO  1590 

1700  IF  G2>J9  THEN  1720 

1710  J9-G2 

1720  NEXT  K1 

1730  J9-J9  MAX  2 HIN  32 

1740  RETURN 
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